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Abstract. We construct projections from HA k (£l), the space of differential 
k forms on f! which belong to L 2 (f2) and whose exterior derivative also be- 
longs to L 2 (Q), to finite dimensional subspaces of HA k (Q) consisting of piece- 
wise polynomial differential forms defined on a simplicial mesh of fi. Thus, 
their definition requires less smoothness than assumed for the definition of the 
canonical interpolants based on the degrees of freedom. Moreover, these pro- 
jections have the properties that they commute with the exterior derivative 
and are bounded in the HA k (Q) norm independent of the mesh size h. Un- 
like some other recent work in this direction, the projections are also locally 
defined in the sense that they are defined by local operators on overlapping 
macroelements, in the spirit of the Clement interpolant. 



1. Introduction 

Projection operators which commute with the governing differential operators 
are key tools for the stability analysis of finite element methods associated to a 
differential complex. In fact, such projections have been a central feature of the 
analysis of mixed finite element methods since the beginning of such analysis, cf. 
However, a key difficulty is that, for most of the standard finite element spaces, 
the canonical projection operators defined from the degrees of freedom are not well- 
defined on the appropriate function spaces. This is the case for the Lagrange finite 
elements, considered as a subspace of the Sobolev space if 1 , and for the Raviart- 
Thomas [18] , Brezzi-Douglas-Marini [7] , and Nedelec [Ml [17] finite element spaces 
considered as subspaces of if (div) or if (curl). For example, the classical continuous 
piecewise linear interpolant, based on the values at the vertices of the mesh, is not 
defined for functions in if 1 in dimensions higher than one. Therefore, even if the 
canonical projections commute with the governing differential operators on smooth 
functions, these operators cannot be directly used in a stability argument for the 
associated finite element method due to the lack of boundedness of the projections 
in the proper operator norms. In addition to the canonical projection operators, 
it is worth mentioning another family of projection operators that commute with 
the exterior derivative. This approach, usually referred to as projection based 
interpolation, is detailed in the work of Demkowicz and collaborators (cf. [8] , [12] , 
[T3] . |14j . [T5]). The main motivation for the construction of these operators was the 
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analysis of the so called p-version of the finite clement method, i.e., the focus is on 
the dependence of the polynomial degree of the finite element spaces. However, as 
in the case of the canonical projection operators, the definition of these operators 
requires some additional smoothness of the underlying functions, so again they 
cannot be used directly in the standard stability arguments. On the other hand, 
the classical Clement interpolant [TT] is a local operator, and it is well-defined 
for functions in L 2 . However, the Clement interpolant is not a projection, and the 
obvious extensions of the Clement operator to higher order finite element differential 
forms (cf. [U [5]) do not commute with the exterior derivative. Therefore, these 
operators are not directly suitable for a stability analysis. 

Bounded commuting projections have been constructed in previous work. The 
first such construction was given by Schobcrl in [19]. The idea is to compose a 
smoothing operator and the unbounded canonical projection to obtain a bounded 
operator which maps the proper function space into the finite element space. In 
order to obtain a projection, one composes the resulting operator with the inverse 
of this operator restricted to the finite element space. In [19] . a perturbation of 
the finite element space itself was used to construct the proper smoother. In a re- 
lated paper, Christiansen [9] proposed to use a more standard smoothing operator 
defined by a mollifier function. Using this idea, variants of Schoberl's construction 
are analyzed in p] Section 5], [3] Section 5], and [10]. The constructed projections 
commute with the exterior derivative and they are bounded in L 2 . Therefore, they 
can be used to establish stability of finite element methods. However, these projec- 
tions lack another key property of the canonical projections; they are not locally 
defined. In fact, up to now it has been an open question if it is possible to construct 
bounded and commuting projections which are locally defined. The projections de- 
fined in this paper have all these properties. The construction presented below 
resembles the construction of the Clement operator in the sense that it is based on 
local operators on overlapping macroelements. 

We will adopt the language of finite clement exterior calculus as in [T] [3] . The 
theory presented in these papers may be described as follows. Let £1 C R™ be a 
bounded polyhedral domain, and let HA k (n) be the space of differential k forms 
u on fl, which is in L 2 , and where its exterior derivative, du = d k u, is also in L 2 . 
This space is a Hilbert space. The L 2 version of the dc Rham complex then takes 
the form 

HA°(Q) 4 HA\n) 4 • • • A HA n (Q). 

The basic construction in finite element exterior calculus is of a corresponding 
subcomplex 

A" 4 A£ 4- -4 AS, 

where the spaces A^ are finite dimensional subspaces of H A k (£l) consisting of piece- 
wise polynomial differential forms with respect to a partition, Tk, of the domain Q. 
In the theoretical analysis of the stability of numerical methods constructed from 
this discrete complex, bounded projections : HA h (fl) —> A^ are utilized, such 
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that the diagram 

7fA°(0) -A HA^tl) -A ••• #A"(fi) 
1"* 1'* 1* 

A d , A 1 d v jj v A n 

^h ' ^ft ' ' ' ' ' ly h 

commutes. Such commuting projections are referred to as cochain projections. The 
importance of bounded cochain projections is immediately seen from the analysis 
of the mixed finite element approximation of the associated Hodge Laplacian. In 
fact, it follows from the results of [H Section 3.3] that the existence of bounded 
cochain projections is equivalent to stability of the associated finite clement method. 
Furthermore, if these projections are local, like the ones we construct here, then 
improved properties with respect to error estimates and adaptivity may be obtained. 

For a general reference to finite clement exterior calculus, we refer to the survey 
papers p] |3] , and references given therein. As is shown there, the spaces K\ are 
taken from two main families. Either A^ is of the form V r A k (Th), consisting of 
all elements of HA k (Vi) which restrict to polynomial /c-forms of degree at most r 
on each simplex T in the partition 7/j, or A k = V~ A k (Th), which is a space which 
sits between V r A k (Th) and V r -iA k (Th) (the exact definition will be recalled below). 
These spaces are generalizations of the Raviart-Thomas and Brezzi-Douglas-Marini 
spaces, used to discretize £T(div) and if (rot) in two space dimensions, and the 
Nedelec edge and face spaces of the first and second kind, used to discretize H(cm\) 
and H(div) in three space dimensions. 

A main feature of the construction of the projections given below is that they 
are based on a direct sum geometrical decomposition of the finite element space. 
In the general case of finite element differential forms, such a decomposition was 
constructed in [2J. However, this is a standard concept in the case of Lagrange finite 
elements. Let Th be a simplicial triangulation of a polyhedral domain 51 G R'\ If 
T is a simplex we let A(T) be the set of all subsimplexes of T, and by A m (T) all 
subsimplcxcs of dimension m. So if T is a tetrahedron in R 3 , then A m (T) are the 
set of vertices, edges, and faces of T for m = 0, 1, 2, respectively. We further denote 
by A(7h) the set of all subsimplices of all dimensions of the triangulation Th, and 
correspondingly by A m (Th) the set of all subsimplices of dimension m. The desired 
geometric decomposition of the spaces V r A k (Th) and V~A k (Th) is based on the 
property that the elements of these spaces are uniquely determined by their trace, 
try, for all / of A (Th) with dimension greater or equal to k. The decompositions 
of the spaces V r A k (Th) established in [2J, is then of the form 

(1.1) PrA k (%)= E k Lr (V r A k (f)). 

/EA(T k ) 

dim f>k 

Here V r A k (f)) is the subspace of V r A k (f)) consisting of elements with vanishing 
trace on the boundary of /. The operator E k r : V r A k (f) — > P r A k (Th) is an exten- 
sion operator in the sense that tr/ °E k r is the identity operator on V r A k (f)). Fur- 
thermore, E k r is local in the sense that the support of functions in E k r (V r A k (f)) 
is restricted to the union of the elements of Th which have / as a subsimplex. A 
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completely analog decomposition 

(1-2) P-A k (%)= E k f -(V-A k (f)) 

dim />& 

exists for the space P~ A k (Th)- 

We will utilize modifications of the decompositions and (|f .2[) to construct 
local bounded cochain projections onto the finite element spaces P r A k (Th) and 
P~ A k {Th)- In the spirit of the Clement operator we will use local projections to 
define the operators tr/ oir k for each / G A(7h) with dimension greater or equal to 
k. To make sure that the projections ir k commute with the exterior derivative we 
will use a local Hodge Laplace problem to define the local projections, while the 
extension operators will be of the form of harmonic extension operators. 

This paper is organized as follows. In Scction[2]wc introduce some basic notation, 
and we show how to construct the new projection in the case of scalar valued 
functions, or zero forms. We also review some basic results on differential forms 
and their finite element approximations. A key step of the theory below is to 
construct a special projection into the space of Whitney forms [20J . i.e., the space 
Vi A k {Th)- In fact, in the present setting the construction in this lowest order 
case is in some sense the most difficult part of the theory, since here we need to 
relate local operators defined on different subdomains. To achieve this we utilize a 
structure which resembles the Cech-de Rham double complex, cf. [4]. In addition 
to being a projection onto the Whitney forms, the special projection constructed 
in Section [3] will also satisfy a mean value property with respect to higher order 
finite element spaces, cf. equation (|3.1[) below. The general construction of the 
cochain projections, covering all spaces of the form P r A k (Th) or P~A k (Th), is then 
performed in Section [4] Finally, in Section [5] we derive precise local bounds for the 
constructed projections. 

2. Notation and preliminaries 

We will use (•, •) to denote L 2 inner products on the domain ft. For subdomains 
D C ri we will use a subscript to indicate the domain, i.e., we write (u,v)d to 
denote L 2 inner product on the domain D. 

We will assume that {Th} is a family of simplicial triangulations of O £ E™, 
indexed by the mesh parameter h = maxTeTh h>Ti where is the diameter of 
T. In fact, hf will be used to denote the diameter of any / 6 A(Th)- We will 
assume throughout that the triangulation is shape regular, i.e. the ratio h^/\T\ 
is uniformly bounded for all the simplices T £ Th and all triangulations of the 
family. Here \T\ denotes the volume of T. Note that it is a simple consequence of 
shape regularity that the ratio hr/hf, for / e A (T) with dim / > 1 is also uniformly 
bounded. We will use [xq, x\, . . . Xk] to denote the convex combination of the points 
x ,Xi,...,Xk e ^- Hence, any / e A k (%) is of the form / = [x , x%, . . . x k ], 
where Xq, x\, . . . , x k 6 Ao(7h). Furthermore, the order of the points Xj reflects the 
orientation of the manifold /. We will let fj £ A k -i{Th) denote the subcomplex of 
/ obtained by deleting the vertex Xj, i.e., fj ~ [xq, . . . , Xj-i, Xj,Xj+i, ■ ■ ■ x k ]- Here 
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the symbol ^ over a term means that the term is omitted. Hence, if j is even then 
fj has the orientation induced from /, while if the orientation is reversed if j is 
odd. 

For each / G A(Th), we let f2/ be the associated macroelement consisting of the 
union of the elements of Th containing /, i.e., 

n f = \J{T\TeT h ,fe A(T)}. 
Fig. 1: Vertex macroelement, n = 2. Edge macroelement, n = 2. 



In addition to macroelements f2/ we will also find it convenient to introduce the 
notion of an extended macroelement f2^ defined for / € A (7/j) by 

«/ = U n v 

SGA (/) 

Fig. 3: The extended macroelement f2^ corresponding to the union of the two 
macroelements tt go (outlined by the thick lines) and f2 ffl , n = 2. 




In the special case that dim / = 0, i.e., / is a vertex, then = ftf. In general, 
if f,geA(%) with. 9 GA(/) then 

flf C fi ff and fi^ C fi}. 

We shall assume throughout that all the macroelements of the form f2/ and f2^, 
for / G A(7/t), are contractive. We let Tf,h denote the restriction of Th to flf, 
while Tf h is the corresponding restriction of Th to Vfj. It is straightforward to 
check that a consequence of the shape regularity of the family {Th} is that the 
ratio |fi^|/|f2/| is uniformly bounded. Furthermore, the coverings {fi/}/6A(Th) an d 
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{fij} /eA(Th) °f the domain both have the bounded overlap property, i.e., the 
sum of the characteristic functions is bounded uniformly in h. 

2.1. Construction of the projection for scalar valued functions. To mo- 
tivate the construction for the general case of k forms given below, we will first 
give an outline of how the projection is constructed for zero forms, i.e. for scalar 
valued functions. The projection 7r° will map the space = HA°(Q) into 

V r AP(7h), the space of continuous piecewise polynomials of degree r with respect 
to the partition Th- The space P r A°(Tf.h) is the restriction of the space V r A a (Th) 
t° Tf,h, an d P r A°(Tf,h) is the subspace of P r A°(Tf\h) of functions which vanish on 
the boundary of fi/, dttf. Of course, by the zero extension the space P r A°(Tf,h) 
can also be considered MS (1 subspace of P r A°(Th)- 

A key tool for the construction is the local projection P° : 7J 1 — > P r (Tf } h), 
associated to each / G A(Th)- If dim / = 0, such that / is a vertex, we define P® 
by P®u G V r A°(Tf,h) as the H 1 projection of u, i.e., P®u is the solution of: 

(P/\l>% = <u,l>n„ 

(dPfU,dv) nj = (du,dv)a J , v G V r (Tf,h)- 

Of course, for zero forms the exterior derivative, d, can be identified with the 
ordinary gradient operator. When 1 < dim / < n, we first define the space 

V r A°(Tf, h ) = {ue V r A°(Tf,h) \tifuG Pr{f) }• 

We then define P®u G P r A°(Tf,h) as the solution of 

{dPju,dv)n f = {du,dv)n f , v G P r A (T f , h ). 

The projection 7r° will be defined recursively with respect to the dimensions of the 
subsimplices of the triangulation Th- More precisely, we will utilize a sequence of 
local operators {7r° n h Y, l n =oi and define = 7r° h . Dropping the dependence on h, 
the operators 7r^ are defined recursively by 

(2.1) ■K° mU = TT° m _ l u+ J2 E °f^f P f( u ~<-i u )^ l<m<n. 

feA m (T h ) 

Here E® : P r {f) — > PrA°(Tf : h) C P r A°(Th) is the harmonic extension operator 
determined by 

(dE° f <j>,dv}n f -0, v eP r A G {Tf, h )M f v = 0, 

and that try £79 is the identity on V r (f). We observe that the dependency of the 
operator E® on the degree r is suppressed. It is a key property that tr g Ej<fr = Q 
for all g G A(Th), dim g < dim/, and g ^ f ■ For the vertex degrees of freedom we 
will use an alternative extension operator. We simply define 7Tq h = ttq by 

n° u = £ £°t r/ P°u = £ ^/°(^)(/) 

f£A (T h ) feA„(T h ) 

where, for any a 6 1, f J a is the piecewise linear function with value a at the vertex 
/ and value zero at all other vertices. Hence, for / G Ao(7fc) we have £2 = E® 
if r = 1. The reason for choosing the special low order extension operator for 
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vertices is not essential at this point, but will be needed later to make sure that the 
projections ~k\ commute with the exterior derivative. 

The key result for the construction above is that the operator 7r° is a projection. 
Lemma 2.1. The operator is a projection onto V r A°(Th)- 

Proof. To see that ir° = 7r° is a projection, we only need to check that if u 6 
V r A°(Th), then for all / 6 A(7h,), tr/7r°it = tr/ti. We do this by induction on 
to, where to corresponds to the dimension of the face / G A(7h). We assume 
throughout that u 6 V r A°(Th)- We will show that the operator has the property 
that 

(2.2) tr/7T^u = tr / u if / G A (7ft) withdim/<TO, 

and since 7r° = 7r° this will establish the desired result. If / G Aq(7/i) then P®u = 
u\n f - By construction, it therefore follows that (|2.2j) holds for m = 0. Assume next 
that (|2.2|) holds for m — 1, where 1 < to < n. It follows that for any / G A m (Th) we 
have tr/(u — 7r° n _ 1 u) G P r {f), and therefore Pj? (u — Tr^^u) = it — Tr.^^jit. It follows 
by construction that tr g 7t^m = tv g tt^ 1 _ 1 u = tv g u for g G A(7/j), with dim g < to, 
while for / G A m (7^) we have 

tvfK^u = trf^^^u + Pj{u - n m _ 1 u) — trf u. 

Therefore, (|2.2[) holds for m and the proof is completed. □ 

It follows from the construction above that the operator 7r° is local. For example, 
for any T G Th we have that (ir® h u)T depends only on u restricted to the extended 
macroelement Define D m T C f2 by 

(2.3) D m , T = U{ £> m _i, T ' I T' G Tf,h, f e A TO (T) }, £> ,t = Of,. 

It follows from (|2.1[1 that (7r m ./iu)|t depends only on w|_o m T . In particular, {^u)\t 
depends only on u\d t , where Dt = P> ni T- 

The operator 7r° satisfies the following local estimate. 
Theorem 2.2. Let T G 7^. TT&e operator 7r° satisfies the bounds 
lk°u|| L 2 (T) < C(||u||£2 (£)t) + /i T ||rfw||z,2(£i r )) 

and 

ll d7r ° M IU 2 (T) < CUduH^^), 
where the constant C is independent of h and T G Th- 
in fact, this result is just a special case of Theorem l5.2l below. so we omit the proof 
here. Of course, due to the bounded overlap property of the covering {DT}reT h 
of fi, derived from the corresponding property of {f^}, global estimates follow 
directly from the local estimates above. 



8 



RICHARD S. FALK AND RAGNAR WINTHER 



2.2. Differential forms and finite element spaces. We will basically adopt 
the notation from [3]. The spaces r P r A k (Th) C HA k (fl) can be characterized as the 
space of piecewise polynomial k forms u of degree less than or equal to r, such that 
the trace, trf u, is continuous for all / £ A(7h), with dim / > k, where we recall that 
the trace, ir/, of a differential form is defined by restricting to / and applying the 
form only to tangent vectors. The space V~A k (Th) C HA k (£l) is defined similarly, 
but on each element T £ Th, u is restricted to be in V~ A k C V r A k . Here, the 
polynomial class V~ A k consists of all elements u of T r A k such that u contracted 
with the position vector x, u_ix, is in T > r A k ~ 1 . Hence, for each k we have a sequence 
of nested spaces 

V{A k (%) c ViA k {T h ) c VzA k {T h ) c . . .HA k (n). 

In particular, p-A°(T h ) = V r A°(T h ), and p-A n (T h ) = V r -iA n (T h ). 

Instead of distinguishing the theory for the spaces V~A k (Th) and V r A k (Th) we 
will use the simplified notation VA k (Th) to denote cither a space of the family 
V~A k (Th) or V r A k (Th)- More precisely, we assume that we are given a sequence 
of spaces VA k (Th), for k = 0,1,..., n, such that the corresponding polynomial 
sequence (VA,d), given by 

(2.4) R^VA°(R n ) -A VA\W n ) -A • • • -A PA n (R n ) -> 

is an exact complex (cf. Section 5.1.4 of [3]). In particular, this allows for combi- 
nations of spaces taken from the two families V~A k (Th) and V r A k (Th)- For any 
/ £ A(T), with dim/ > k, the space VA k (f) = ti f VA k (T h ): while VA k (f) = 
{v £ VA k (f) | tigfv = 0}. The corresponding polynomial complexes of the form 
(VA(f),d) are all exact. Furthermore, the complexes with homogeneous boundary 
conditions, (VA(f),d), given by 

(2.5) VA°(f) PA^f) -i> ... -1> ■pA dbn f (/) -)■ M 
are also exact. 

We recall that the spaces VA k (7~h) admit degrees of freedom of the form 

(2.6) [ttfuAr,, r)€V'(f,k), f £ A(T h ), 
Jf 

where V(f,k) C A dlm *~ k (f) is a polynomial space of differential forms and the 
symbol A is used to denote the exterior product. These degrees of freedom uniquely 
determine an element in VA k (Th), (cf. Theorem 5.5 of [3]). In fact, if 

VA k (%) = V~A k (Th), then V'(f, k) = P r+fc _ dim/ _ 1 A dim/ - fc (/), 

while if 

TA k (%) = V^iTn), then P'(f, k) = ^+ fe _ dim/ A dim ^ fc (/). 

If v £ VA k (f), then v is uniquely determined by the functionals derived from 
V(f,k). Furthermore, any v £ VA (/) is uniquely determined by V'(g,k) for all 
g £ A(/). In particular, if dim f < k then V'(f, k) is empty, while V'(f, k) is always 
nonempty if dim / = k. For dim/ > k the set V(f,k) can also be empty if the 
polynomial degree r is sufficiently low. 
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The local spaces VA k (Tf,h) and VA k (Tj h ) are defined by restricting the space 
VA k (Th) to the macroelements f2/ or SVj. It follows from the assumption that fi/ 
and £Yj are contractive, that all the local complexes (VA(Tfji),d) and (VA(Tf h ), d) 

are exact. The same holds for the subcomplexes (VA(Tf,h),d) and (P A(Tf h ) 7 d) 7 
corresponding to the subspaces of functions with zero trace on the boundary of the 
macroelements . 



For a given triangulation Th, the spaces of lowest order polynomial degree, 
Vi A k (Th), i-e., the space of Whitney forms, will play a special role in our construc- 
tion. The dimension of this space is equal to the number of elements in A k (Th), 
and the properties of these spaces will in some sense reflect the properties of the 
triangulation. Therefore, this space will be used to transfer information between 
different macroelements, cf. Section[3]below. For k = this space is just V\A k (Th), 
the space of continuous piecewise linear functions. The natural basis for this space 
is the set of generalized barycentric coordinates, defined to be one at one vertex, 
and zero at all other vertices. It follows from the discussion above that the degrees 
of freedom for the space Vy A k (Th), < k < n, are f.u for all / G A/ £ (7/ l ). In 
fact, if / = [xq,x\, . . .Xk] G A.fc(77i), wc define the Whitney form associated to /, 
4>) G VyA k {T h ), by 

fc 

4> k = ^2(-iy\id\ A • • • A d\i A • • • A d\ k , 

i=0 

where Ao, Ai, . . . , Afc are the barycentric coordinates associated to the vertices x%. 
The basis function <p k reduces to a constant k form on /, i.e., tr/ <p k G PoA k (f), and 
it has the property that tr s cj) k = for g G Afc(7fe), g ^ /• In fact, if vol/ G VoA k (f) 
is the volume form on /, scaled such that vol/ = 1, then 

tr/ 4>) = (k\)-\o\ f: 

cf. [1] Section 4.1]. Furthermore, the map vol/ — > £ k vo\f — k\<j)h defines an 
extension operator £ k : P A k (f) -> V 1 ~A k (Tf,h) for any / G A k (Th)- We observe 
that the operators £ k are natural generalizations of the piecewise linear extension 
operators £/, introduced above for scalar valued functions. In fact, any element u 
of A k (Th) admits the representation 

(2.7) u= Y, (f tr f u 



We finally note that it follows from Stokes' theorem that if / = [xq, x\, . . . , x k +i], 
and it is a sufficiently smooth k form on /, then 

/fc-f-i /> 
du = Y(- 1 ) j l ^«» 
j—0 fj 

where fj = [xo, . . . , Xj—i, Xj, Xj+\, ■ • ■ Xk-\-i\- Here the factor (— \) 3 enters as a 
consequence of orientation. 
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3. A SPECIAL PROJECTION ONTO THE WHITNEY FORMS 



Recall that the purpose of this paper is to construct local cochain projections 
7r£ which map HA k (fl) boundedly onto the piecewise polynomial space VA k (Th)- 
Furthermore, in the construction of 7r° given above, the construction of tr/o7r° is 
based on a local projection, Pj?, defined with respect to the associated macroelement 
0/. Therefore one might hope that all the projections ir^ have the property that 
tr/ o7r^ is defined from a local projection operator defined on fi/ for / £ A(7h), 
dim/ > k. However, a simple computation in two space dimensions, and with 
VA k (Th) = Vy A k (Th), will convince the reader that if / = [#0,2:1] £ Ai(7fc), then 

f Xl d 



f Jx ds 

and the right hand side here clearly depends on u restricted to the union of the 
macroelements associated to the vertices xq and xy. Therefore, JjtrfKydu = 
ff tr/ dir^u must also depend on u restricted to the union of these macroelements, 
and this domain is exactly equal to the extended macroelement fi/. This motivates 
why the extended macroelements, f2j, for / £ Afc(7ft,), will appear in the construc- 
tion below. In fact, a special projection operator, i?^ : HA k (Vl) — > Vy A k (Th) C 
VA k (Tii), will be utilized in the construction of ir^ to make sure that 

tx f % k h du = tr f dir k - 1 u= / traf^u, 

Jf Jf Jdf 



for all / £ A k (%). 



The operator will commute with the exterior derivative, and it is a projec- 
tion onto Vi A k (Th)- Therefore, in the case of lowest polynomial degree, when 
VA k (Th) — Vy A k (Th), we will take tt^ = R 1 ^. However, another key property of 
the operator R k is that in the general case, when Vy A k {Th) is only contained in 
VA k (%), we will have 

(3.1) tr/ R k h u = tr/ u, f £ A k (%), u £ TA k (T h ), 

i.e., the operator R k preserves the mean values of the traces of function in VA k (Th) 
on subsimplexes / of dimension k. The rest of this section is devoted to the con- 
struction of the operator R k , and the derivation of the key properties given in 
Theorem 13.61 below. 



3.1. Tools for the construction. Throughout this section the dependence on 
the mesh parameter h is suppressed in order to simplify the notation. To define 
the special projection R k onto the Whitney forms, VyA k (T), we will use local 
projections, Q k , defined with respect to the extended macroelements il/. We define 
the projection Q k f : HA k (tt e f ) —} VA k {Tf) by the system 

(Q k f u,dT) a . = (u,dr)n } , r £ VA k ~ x {Tf\ 
(dQ k f u,dv)^ f = (du,do)m, v £ VA k (T f e ). 
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For k = 0, the first equation should be replaced by a mean value condition, so that 
Q°j: = P®. This system has a unique solution due to the exactness of the complex 
(V A(7J) , d) . Furthermore, by construction we have 

(3.2) Q)du = dQ k f x u, 0<k<n. 

We will also find it useful to introduce the operator Q k f _ : HA k {VL e f ) — » VA k ~ l {Tf) 
defined by the corresponding reduced system 

{Ql_u, dr) n} = 0, r G VA k - 2 (T f %), 

(dQ k ft _u,dv) U} = <u,dw)nj, v € VA 1 " 1 ^). 

As a consequence, the projection Q k can be expressed as 

(3.3) Q k f = dQ k ft _ + Q)+}d. 

To make this relation true also in the case when k = and / <E Ao(T), the operator 
dQ°f _ should have the interpretation that dQ® _u is the constant / n u A voln^ on 
Qf, where voIq / is the volume form on f2, restricted to O j and scaled such that 
In, vo1 ^ = I- 

To motivate the rest of the tools we need to for our construction, consider again 
diT u in the special case when VA k {T) = P-fA^T). To obtain a commuting 
relation of the form dir°u = w du, we have to be able to express dir°u in terms 
of du. However, using the notation just introduced, we have 

dn a u= ^2 \( u Avolfig) +ti g (Qg_du) dE g vo\ g . 

The second part of this sum is already expressed in terms of du. By combining the 
contributions from neighbouring macroelements we wil see that also the first part 
of the right hand side can be expressed in term of du. If / = [xq, x\\ G Ai(7~), wc 
have 

/ tr / J2 (/ uAvoln a )d£ fl vol fl = / u A (vob sl - voln ao ), 

f ff6A (T) ■« n "f 

where gi = [xi]. Furthermore, voln 91 — V0 '^ 90 e PoA n (Tj) = V^A n (Tf), and with 
vanishing integral. As a consequence, there exists zj G 'P 1 "A" _1 (7j : ) such that 
dzj = volo so — vo 'f2 gi j an d by integration by parts 

/ tr/ ( u A vo\ci2)d£g\/o\g = — / u A dz j = / duAzj. 

f gGAo(T) f / / 

By utilizing the representation (|2.7|) . wc therefore obtain 



E 



/ u Avo\n g ^jd£ g vo\ g = ^ ^ / du A zjjE k \/o\f. 



9 eA (T) /eAi(Th) Jo / 

This discussion shows that to construct local cochain projections, we must utilize 
relations between local operators defined on different macroelements. To derive the 
proper relations, we introduce an operator 

5: V-A k (T g e )^ V^A k (T f e ). 

ff eA m (T) /eA m+1 (T) 
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Xq , . . . , Xm-\-l 

] € A m +i(T), then the component (5u)f of is defined by 

m+l 

(H/ = E(- 1 )^> 

3=0 

where, as above, /j = [xo, . . . , Xj-i, ij, x^+i, ■ • ■ x m +i], and it/^. the corresponding 
component of u. We will also consider the exterior derivative d as an operator 
mapping © s6Am(r) V^A k (T g e ) to © seAm(r) A^ 1 (T*) by applying it to each 
component. Hence, the two operators doS and dod both map © 9e A m (T) ^ k 0~g) 
into ©/ e A m i(T) T^i ^ h+1 (Tf) ■ In fact, we have the stucture of a double complex 
which resembles the well-known Cech-de Rham complex, cf. [1]. The following 
two properties of the operator S are crucial. 

Lemma 3.1. 

d o S = S o d, and 5 o 5 = 0. 



Proof. It follows directly from the definition of 6 that for / = [xo, . . . , x m +i] G 
A m+ i(T) 

m+l 

(do5u) f = (5odu) f = ^(-l) J cZw /3 . 

3=0 

If we further denote by the subsimplex of / obtained by deleting both Xi and 
Xj, then 

m+l 
3=0 

m+l j — 1 m+l 



^(-i)^(-i)^- 53 (-i)^ 



3=0 



i=0 



0. 



since for each i,j = 0, ...m + 1, with i ^ j, the term u/ y appears exactly twice 
with opposite signs. □ 



The construction of the projection R k = R k will depend on local weight func- 
tions, z) e P 1 "A n - fe (7} e ) for / G A fc (r). In particular, the function z° G V A n (T f ) 
for / G Ao(T) will be given by z'j = voln / - For k = 1,2, ...,n, the functions 
2:^ G Vi A n ~ k (Tj) are defined recursively to satisfy the conditions 

(3.4) <fe* = (-l)*^*- 1 )/, 
and 

(3.5) (4, dr) n? =0, r G PfA"-^ 1 ^/), 

for any / G Afe(T). We will not give an explicit construction of the functions zH. 
However, we have the following basic result. 

Lemma 3.2. The weight functions z k j G 'P-f A n ~ k (Tj) exist and are uniquely de- 
termined by Zt and the conditions (|3.4[) and (|3.5I) • 
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Proof. We establish the existence of the functions Zj by induction on k. Let / = 
[x ,xi] G Ai(T). Then 

(5z°) f = (z° fl - z%) = vol n/l - voln /0 , 

which implies that J ne (Sz°)f = 0. Hence, by the exactness of the complex 

PfA"- 1 ^/) A V A n {Tf) -> R, 
there exists zj- € A" _1 (7j : ) satisfying dzj = — (5z )/. Next, assume we have con- 
structed z*- 1 g © /e A fc _ 1 m^r An-fc+1 (7f ) such that = (- 1 ) k ~ 1 ( Szk ~ 2 )f 

for all / G Afc_i(T). From Lemma lOl we obtain 

(doS)z k - 1 = {5od)z k - 1 = (-l) k - 1 (6o5)z k - 2 = 0, 

and for each / G A/ £ (7") the complex (d, A(7^ e )) is exact. Therefore, we can 
conclude that there is a z k f G "P 1 ~A™- fe (A fe (T f e )) such that (pTl]) holds. This com- 
pletes the induction argument. Finally, we observe that it is a consequence of (|3.5p 
and the exactness of the complex (d, A(Tf)) that each function z k is uniquely 
determined. □ 



We will use the functions z k to define an operator M k := M k : HA k (il) — > 
V^A k (T) by 

M k u= (/ uAzfj£j'vo\f. 

Note that M k u is a generalization for £;-forms of the expression 

Y (/ uAvol n/ )f/vo]/. 

appearing above in the case of zero-forms. It follows from the construction of the 
functions z k that the operator M k commutes with the exterior derivative. 

Lemma 3.3. For any v G HA k ~ 1 (fl) the identity dM k ~ 1 v = M k dv holds. 



Proof. We have to show that 

E (/ vAz k g -')d£ k -\o\ g = £ (/ dvAz k )£ k vo\ f 

seAi-ifT) Jn s /eA fe (r) JO / 

for any v G HA k ~ 1 (n). Since both sides of this equation are elements of 'P ] ~A fe (T), 
we need only check that the integrals of their traces are the same over each / = 
[xq.Xi, . . . ,Xk] G Afc(T). Now it follows from the properties of the extension oper- 
ators Eh that the integral of the right hand side is simply f Qe dv A z k , while (|2.8|) 

implies that the corresponding integral of the left hand side is 
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where the last identity follows by (|3.4|) . However, by integration by parts (cf. [U 
Section 2.2]), utilizing that trgsij Zf = 0, we have 

/ vAdz) = (-l) fc / dv Az k f , 

and this completes the proof. □ 

We now define an operator S k = S k : HA k (Q) — > P^A k (T) recursively by 
5° = M° and 

S k u = M k u + ( / tr fl[ 7 - S^Ql^d^-^olg, l<k<n. 

We recall that the operator Q k t _ is a local operator with range PA{Tf). However, 
by an inductive argument, it follows that the composition try oS k is a local operator 
mapping HA k (n e f ) into V A h (f). Therefore, the operators S k are indeed well 
defined. 

A key property of the operator S k is the following result. 
Lemma 3.4. For any v G PA k ~ 1 (T) the identity 

^S k dv = J^dv, f€A k (J). 

holds. 

Proof. The proof goes by induction on k. For k = the space dP (T) should be 
interpreted as the space of constants on fi, and since 5° = M° reproduces constants 
the desired identity holds. 

Assume next that k > 1 and that the desired identity holds for k — 1. By 
utilizing the result of Lemma l3~3l we obtain that the commutator S k d — dS k ~ 1 has 
the representation 

S k dv~dS k ~ 1 v = ( f tegil-S^Q^-dvyS^volg, v e HA k ~ 1 (tt). 

We recall that Q k _dv = Q k ~ 1 v—dQ k Z L v, and by the induction hypothesis f tr g (I — 

S^dQ^v = 0. Therefore, for any v e i/A fc 1 (f2) the commutator above can 
be expressed as 

(3.6) S k dv-dS k - 1 v= J2 ( f tt g [I - S^C^-^dS^volg. 

However, since is a projection onto PA k ~ 1 (T g e ). it follows from (|2.7jl that 

dS k ~ 1 v = J2 ( ftv g S k - 1 Q*- 1 v)d£%- 1 vo\ g , BGPA'-^T). 

By restricting to a function u € PA k ~ 1 (T) equation (|3.6[) therefore reduces to 
S*(fo = ( j tx g v\d£ k - l M0\ g . 



LOCAL BOUNDED COCHAIN PROJECTIONS 



15 



By integrating this representation over any / = [xq, x\, . . . , Xk] € Afc(T) we obtain 



k 

f S k dv = J2(-iy f tr fi v= f 

Jf j = Q J fj Jf 



dv, 



j= 

where the final identity follows from (|2.8p This completes the induction argument. 

□ 



As a direct consequence of the proof above, we have. 
Lemma 3.5. The identity (|3.6[) holds for k = 1,2, ... ,n and all v <E HA k ~ 1 (fl). 

3.2. The projection Rfi. For each k, < k < n, the operator R k = i?^ : 
HA k (fl) -)• V7A k {T) is defined by 

R k u = S k u + (/ tr /t J - S k ]Q)u\e k f M0\ } . 

f£A k (T) J f 

Recall that the operator S k is local in the sense that tr/o,S ,fc can be seen as a 
local operator mapping HA k {Vl ( j) onto PA k (f). It is immediate from this and 
the properties of of the projection Q k , that try oR k also is local. In fact, for any 
T e T, {R k u)\r only depends on u\n%.- Furthermore, if / e A (Th) then Q° f = Pj?. 
Therefore, it follows that for k = the operator P° is identical to the operator 
7Tq h , used in the construction of the projection 7r° in Section [2 . 1 1 above . 

The key properties of the operator i?^ are given in the theorem below. 

Theorem 3.6. The operators R k : HA k (£l) — > P-f A k (Th) are cochain projections. 
Furthermore, they satisfy property (|3.1[) . i.e., 

J tv f R k u = J try u, / G A k {%), u G PA fe (7I). 



Proof. As above we use a notation where we suppress the dependence on h. It is 



a consequence of the projection property of the operators QH that if u € PA fe (T) 



then 



R k u= (/ tr / M )^ vo1 / 



/eA fe (T) 

However, this implies the identity (|3.1[) . and an immediate further consequence is 
that R k is a projection onto V7 A k (Th) ■ 

It remains to show that R k commutes with the exterior derivative. From the 
definition of R k and Lemma l3~5l we have 



dR k u = dS k u + { [ : tr /t J ~ S k ]Q k f u\d£)\io\ f = S k+1 du. 

feA k (T) Jf 

However, S k+1 du = R k+1 du since 

R k+l du - S k+1 du = ( / tr /t J - S k+1 ]Q k+1 du)£ k+1 vo\ f 



f£A k+1 (T) 
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= J2 ( / tr /t J - S k+x ]dQ)+}d^£ k f + \o\f = 0, 
/eA fc+1 (T) J f 

where the last identity follows from Lemma f3T4l □ 

The operators R 1 ^ introduced above are local operators in the sense that (R k u)\T 
only depends on u|n=, for any T € Tn- Furthermore, for any fixed h the operator 
is a bounded operator on HA k (Q). The discussion of more precise local bounds 
is delayed until the final section of the paper. 

4. Construction of the Projection: The General Case 

We finally turn to the construction of the projections in the general case, in 
which VA k (Th) denotes any family of spaces of the form V~A k (Th) or V r A k (Th), 
such that the corresponding polynomial sequence (VA k , d) , given by (|2.4[) is an exact 
complex. In particular, the Whitney forms, Vy A k (Th), are a subset of VA k (%), 
and in the special case when VA k (Th) = Vy A k (Th) we will take -k\ to be the 
operator R\ constructed above. 

In the construction we will utilize a decomposition of VA k (Th) of the form 

(4.1) VA k (T h )= £ k (P A k (f)) + E k {VA k (f)), 

feA k (T h ) /EA(T h ) 

dim f>k 

where £ k is the extension operator defined in the previous section, mapping into 
the space of Whitney forms, while E k is an harmonic extension operator mapping 
into VA k (Tf, h )- Furthermore, the space "PA fe (/) = VA k (f) if dim / > k, while 

VA k (f) = {u 6 VA k (f) | J f u = °}, if dim/ = fc. 

The decomposition (|4.ip can be seen a modification of the more standard decom- 
positions and (|1.2p in the sense that we arc utilizing the special extension, £ k , 
for the constant term of the traces on /, when dim / = k. The existence of such a 
decomposition of the space VA k (Tf) is an immediate consequence of the degrees of 
freedom (|2.6|) . 

As in the case k = 0, cf. Section [2TTT the projection 7r£ will be constructed from 
a sequence of operators 7r^ h , where ir k = n k h . The operators 7r^ h are defined by 
a recursion of the form 

(4.2) 7T* >fc = 7T* + Ek f°^S oP f\ I -^m-i,hl k<m<n, 

f£A m (T h ) 

where the operators P k are local projections defined with respect to the macroele- 
ments f2/, generalizing the operators Pj? introduced in Section 12.11 Furthermore, 
the operator tt\_ 1 h will be taken to be the operator R\ defined in Section [3] above. 
Hence, to complete the definition of 7r^, it remains to give precise definitions of the 
local operators Eh and P k . 
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4.1. Extension operators. As above, to simplify the notation, we suppress the 
dependence on h throughout the discussion. The extension operators E k are gener- 
alizations of the harmonic extension operators E® used for zero forms in Section fe.ll 
Let us first assume that / G A(T) such that / is not a subset of the boundary of 
Qf. In this case, the harmonic extension E k maps VA k (f) to VA k (T}), where 
< k < dim/. More specifically, we let E k cf> be characterized by 

\\dEfy\\ LHQf) = mf{\\dv\\ L 2 (Qf) | v e VA k (T f ), tr f v = 4>}. 

We should note that it is a consequence of the degrees of freedom of the spaces 
VA k (Tf) and VA k (Tf) that there are feasible solutions to this optimization problem. 
As a consequence, an optimal solution exists. However, the solution is in general 
not unique. The solution is only determined up to adding functions w in 'PA fc (7/ j / t ) 
satisfying dw = on fij and trfw — 0. Therefore, to obtain a well defined extension 
operator, we need to introduce a corresponding gauge condition. Hence, for any 
4> E VA k (f) we let E k f 4> G VA k (Tf) be the solution of the system 

{E k 4>,dr)n f =0, r e N(tr f ;VA k -\Tf)), 
{dE k 4>,dv) Ut =0, v E N(tv f ;VA k (T f )), 

and such that tvjoE k is the identity on VA k (f). Here iV(tr/;A) denotes the 
kernel of the operator tr/ restricted to the function space X. A key property of the 
extension operators E k is that they commute with the exterior derivative. 

Lemma 4.1. Let f E A(T). The extension operators E k : VA k (f) -> VA k (T f ) 
are well defined by the system (|4.3|) for k = 0, 1, . . . , dim /, and for k > 1 we have 
the identity 

(4.4) E k f d4> = dE^ 1 ^, (/) <E f>A k ~ 1 (f). 

Moreover, the kernel of d restricted to N(tr f ;TA k (Tf)) is dN(tXf]VA k -' L {Tf)). 

Proof. For k = the first equation in the system (|4.3j) should be omitted. The 
kernel of d restricted to N(tif,PA°(Tf)) is just the zero function, and E'jcj) is 
clearly uniquely determined by the second equation and the extension property. 
We proceed by induction on k. 

Assume that the statement of the lemma holds for all levels less than k. We 
first establish the characterization of the kernel of d, restricted to iV(tr/; VA k (7~f))- 
Assume that u £ N (ti f ; V A k (Tf)) satisfies du = 0. Then, by the exactness of the 
complex (VA(Tf),d), u = dr for some r S VA k ~ 1 (Tf). Furthermore, dtrjr — 
tvfdr = tr f u = 0. If k = 1 this implies that t G N(tv f ;VA (Tf)). For k > 1 it 
follows from the exactness of (VA(f),d) that there is a (j> G VA k ~ 2 (f) such that 
d(j) = tvf t. However, the function 

a = t - dE k ~ 2 4> = t - E k - X d4) 

G N{ttf,VA k - 1 (Tf)) and satisfies da = u. Hence the complex (iV(tr/; VA{Tf)), d) 
is exact at level k in the sense that <i/V(tr/; VA k ~ 1 (Tf)) is the kernel of d restricted 
to N(tT f ;VA k (Tf)). 
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Consider a local Hodge Laplace problem of the form: 

(a,r) Qf - (u,dr}n f = 0, r G N(tr f ; PA k ^(T f )), 

(4-5) 

(du.^Oj. + {du,dv) Qf = 0, we N(tr f ;PA k (T f )), 

where the unknown (cr,it) € iV(tr/; PA fc-1 (7/)) x PA k (T f ) 1 and with tr/ u — (f> G 
PA k (f). Since the complex (N(trf,PA(Tf)),d) is exact at level fc, it follows from 
the abstract theory of Hodge Laplace problems, cf. for example [3J Section 3], that 
the system (|4.5[) has a unique solution. Furthermore, by the exactness of the same 
complex at level k — 1, a = 0. Hence, u and E k <f> satisfy the same conditions, and 
the uniqueness of E k cj) follows by the uniqueness of u. 

Finally, to show the identity (|4.4[) we just observe that for any </> e PA k ~ 1 (f), 
the pair (a, u), with a = and u = dE k ~ l (f> G VA(Tf), satisfies the system (|4.5|) 
with tr/ dE k ~ l § = dtrf E^ 1 ^ = dip. By uniqueness of such solutions we conclude 

that dEf^ 1 ^ = E k d(f>. This completes the induction argument, and the proof of 
the lemma. □ 



If g G A(7/), with k < dimg < dim / and g ^ /, then tr g £j?(/> = 0. In the case 
that / C <9f2, we will also have that / C d£lf . In this case, the definition of the 
operator E k should be properly modified, such that E k <fi is not required to be in 

"PA fe (7/), but only required to be zero on the interior part of dtlf. The key desired 
property is that the extension of E k (f> from £lf to f2, by zero outside flf, is in the 
global space PA k {T h )- 

It is a consequence of the decomposition (|4. 1[) that any element u of VA k {T) 
is uniquely determined by its trace on /, tr/u, for all / G A (7~) with dim / > fc. 
Furthermore, if u is an element of the subspace given by 

(4.6) £ k (P A k (f)) + E k (PA k (fj), 

feA k (T) /SA(T) 

A:<dim /<m 

then u is determined by tr^-u for all / G A(7~) with fc < dim/ < m. A key 
observation is the following. 

Lemma 4.2. Assume that u G PA k {T) belongs to the subspace given by (|4.6[) . 
where k < m < n. Then its exterior derivative, du, belongs to the corresponding 
space 

£ k+1 (P A k+1 (f)) + E k+l {PA k +\f)). 

/GA fc + 1 (T) /EA(T) 

&+l<dim /<m 

Proof. It follows from the fact that (P^A(T),d) is a complex that d£ k \io\ g G 
©/eA fc+1 (r) f / +1 ( 7, ° Afe+1 (/)) for an y -9 e M 7 ")' Furthermore, if .g G A(T) and 
dimg > fc then (|4.4j) implies that = E k+1 d(j) for any G PA k (g). As a 

consequence, it only remains to check terms of the form dE k cj>, where <j> G PA k (g) 
and dim g = k. 
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Note that dE k (f> is identically zero outside f2 9 . Furthermore, consider any / 6 
A fc+ i(T), with g G A fe (/). Then fl f C fl g and the space N(tv f ;'PA k (Tf)) can 
be identified with a subspace of N(tr g ;VA k (T g ))- Therefore, it follows from the 
definition of E k cj) that 

{dE k g <j>,dv) n , = 0, v G N(tr f ;VA k (T f )), f G A fc+1 (T), .9 G A fc (/). 
However, this implies that 

dEfre E k +\VA k+1 {f)) 

/eA fc+l( r) 

geA fe (/) 

= E k+ \V Q K k+1 {f))+ 4 +1 (PA fc+1 (/)). 

/eA fe+1 (T) /GA fc + 1 (T) 

geA k (f) 9 eA fe (/) 
This completes the proof. □ 



The harmonic extension operator just discussed is the one we will use in the 
construction of the local cochain projection 7r fc , cf. (|4.2[) . However, in the theory 
below we will also utilize an alternative local extension, defined with respect to 
spaces VA k (T f ) instead of VA k (Tf). For < k < n the operator E k : VA k (f) 
VA k (T) is defined by the conditions 

{E k 4>,dr)n s =0, r e N(tr f ;VA k -\T f )), 

{dE k cp,dv)n s =0, V G^(tr /; ^A fc (r/)), 

in addition to the extension property try oE k <fr = (j> for all (j> G VA k (f). In complete 

analogy with the discussion for the operators E k above, by utilizing the exactness 

of the complex (VA(Tf),d) instead of the exactness of (VA(Tf),d), we can conclude 
with the following analog of Lemma 14.11 

Lemma 4.3. Let f G A(T). The extension operators E) : VA k (f) ->■ VA k (T f ) 
are well defined by the system ()4.7j) for k = 0, 1, . . . , dim /, and for k > 1 we have 
the identity 

E k f d<p = dE k ~ 1 <p, 4> S VA k ~ 1 (f). 
Moreover, the kernel of d restricted to N(tr f ;PA k (T f )) is dN(ti f ; VA k - 1 (T f )). 

4.2. Local projections. Let / G A(T) and recall the definition of the spaces 
VA k (f) given above, as VA k (f) if k is less than dimension of /, and as the subspace 
of VA k (f) consisting of functions with zero mean value if k = dim /. Hence, as an 
alternative to (|2.5p . we can state that the complex 

0^PA°(/) -A VA 1 {f) pA dim f(f)^0 

is exact. In particular, this means that the first operator, d = d , is one-one and 
the last operator, d = d m -'~ 1 , is onto. In order to define the local projections P k , 

appearing in (|4.2[) . we will use the spaces V(f) to introduce proper local spaces, 
VA k (T f ). For < k < dim / these spaces lie between VA k {T f ) and VA k (T f ), i.e., 

VA k {T f ) G VA k {T f ) G VA k (T f ). 
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More precisely, for < k < dim / the space PA k (Tf) is defined by 

PA k (T f ) = {u £ PA k (T f ) | tv f g PA k {f) }, 

while we let PA k (T f ) = PA k (T f ) for dim / < k < n. We note that for k = this 
definition is consistent with the definition of the space P r A°(7/) used in Section lOl 

We observe that dPA k (T f h ) C PA k+1 (7/ /»)• In other words, (PA k (T f ), d), given 

by 

0^PA°(T/) -A PA\T f ) ••• -A PA"(7»->0, 
is a complex. We also have the following: 

Lemma 4.4. TTie complex (PA k (Tf),d) is exact. 

Proof. Let m = dim/, and assume that u £ PA k {Tf) satisfies du = 0. We need 
to show that there is a a £ PA k ~ 1 (Tf) such that da = u. For k > m + 1 
this follows from the exactness of the complex (PA(Tf),d). Assume next that 
k < m. Since d try u = tifdu = 0, it follows from the exactness of the com- 
plex (PA(f),d) that there is G PA k ~ 1 (f) such that d<j) = tr f u. Therefore 
u - dE)- 1 ^ is in N(tr f ,PA k (T f )) and d(u - dE k f l (t>) = 0. By Lcmma[01 there 
is a r G iV(tr/, PA k ~ 1 (Tf)) such that dr = u - dE k r x §. Hence, the function 
<t = r + satisfies tr/ cr = cp £ PA k - 1 (f). So a £ PA fc " 1 (7/) and dcr = m. 

Finally, we have to consider the case when k — m + 1. The exactness of the 
complex (PA(Tf), d) and the assumption du = implies that there is t G PA m (Tf) 
such that c?r = it. Furthermore, the exactness of (PA(f),d) implies that there is a 
4> £ PA m ~ 1 (f) such that d<j) = tr/ r. The function a = r — dEJ 1 ' 1 ^ has vanishing 
trace on /. Therefore, it is in PA m (Tf), and da ~ u. □ 

We are now ready to define a local projections P k : HA k (ftf) — > PA k (Tf) 
satisfying 

(P}u,dr)a f = {u,dr}n f , r GPA fe - 1 (7», 
(dPfu,dv)n f = (du,dv) nf , v £ PA k (T f ). 

The operator P* is a well defined projection onto PA k (Tf) as a consequence of 
Lemma I4T41 When k = 0, the space d'PA _1 (7/) should be interpreted as the space 
of constants on ffy, such that P° is exactly the projection defined in Section l!Q1 

With this definition it is straightforward to check that the projections P k commute 
with the exterior derivative, i.e., 

(4.8) P k f du = dP^u, 0<k<n. 

4.3. Properties of the Operators 7r^. The definitions of the operators E k and 
P k given above complete the construction of the operators 7r fc = 7r^ given by the 
recursion (|4.2[) . Here we shall derive two key properties of these operators, namely 
that they are projections onto PA k (T) and that they commute with the exterior 
derivative. It is also clear from the construction that the operator is local, 
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and, for each triangulation T — Tin 71 ~h 1S wen defined as an operator on HA k (VL). 
However, the derivation of more precise bounds will be delayed until the next 
section. 

We recall that the recursion (|4.2[) is initialized by choosing 7T*_i = R k , i.e., 
the special projection onto the Whitney forms constructed in Section [3] above. 
Therefore, we obtain from Theorem 13.61 that 

(4.9) dTtl_ 1 u = Tr k+1 du, k = 0, l,...,n-l, 

and for k = the two operators tt^_ 1 and 7Tq are the same. Furthermore, for 
functions in T J A k (T~) the operator n k _ 1 preserves the integral of the trace over all 
subsimplexes of dimension fc, i.e., 

(4.10) Jtr f *£_ 1 u = Jixfu, / G A fe (T), u G VA k {T). 

In other words, if u € VA k (T), then (u - n^u)^ G VA k (T f ) for / G A fc (T) and 
fc > 1. 

We observe that it follows from (|4.2p and the properties of the extension operators 
E k , that if / G A m (%), with m > fc, then 

(4.11) try TT^U = tr/(7T* _!« + P^[u ~ 7T* _ lU ]). 

On the other hand, 

(4.12) tr g 7if n u = tr ff TT^jti, g G A(T), fc < dimg < m. 
These observations are the key tools to show that 7r fc is a projection. 
Theorem 4.5. The operators tt^ are projections onto VA k (Th)- 

Proof. Assume throughout that u G VA (T). We have to show that Tr k u = u. We 
will argue that 

(4.13) tryTT^u = tr/u, if / G A(T), fc<dim/<m, 

for m = fc,fc + l,...,n. This will imply the desired result, since functions in 
VA k (T) are uniquely determined by their traces on / G A (T). We will prove 
(|4.13|) by induction on to. Recall that it - G VA k (Tf) for any / G A fc (T). 

As a consequence, (w — 7r£ _ : u) = M — Tr^'^-it, and therefore (|4. 13[) . with to, = fc, 
follows from (|4. 1 1|) . 

Next, if (|4.13p holds for m replaced by to — 1, then (|4.12p implies that tr g TT k n ii = 
tr g tt!^ 1 _ 1 u = tr g u for all g G A(T), with fc < dim/ < m. So it only remains to 
show the identity (|4.13j) for / G A m (T). However, for each / G A, n (T), we have 
(it - Km-iu)\n f G VA k {Tj). Hence Py (u _ ^m-i") = ( u _ ^m-i^ln^ and then 
(|4.11[) implies that tifir k n u = tryu. We have therefore verified that the operator 
7r^ satisfies property (|4.13j) . This completes the proof. □ 



To show that the projections tt^ are cochain projections, the following observa- 
tion is useful. 
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Lemma 4.6. Assume that < k < n and that u e HA k For any f 6 Afc(T) 

the function d(7T k Zlu — ^Z^OIn/ £ VA k (7}). 

Proof. The function e = d(7r^Ju — 7r£~2 u ) is obviously in VA k (Tf). Therefore, it 
only remains to show that try e = 0. If / = [xq, x\, ■ ■ ■ , then it follows from 
the definition of n^zl and (|2T5|) that 

/ tr/ c = / tx h P /~> - = °' 

/ j — Q fj 

Here the last identity follows since for dim/j = fc — 1, the projection Py. _1 projects 
into a space of functions of mean value zero on fj . □ 

We conclude with the final result of this section. 

Theorem 4.7. The operators 7r£ are cochain projections, i.e., dirZ~ 1 = n^d for 
k = 1, 2, . . . , n. 

Proof. We will prove that for u £ HA (f2), and 1 < k < n, 

(4.14) tT f TT^ n du = tTfdTT k r Z 1 u, if / e A(%) , k < dim / < m, 

for m = fc, fc + 1, . . . , n. As above, the case to = n implies the desired result. We 
note that it follows from Lemma H. 2 1 that if (|4.14[) holds for any fc < m < n, then 
ir k n du = dir 1 ^Z 1 u. 

The identity (|4.14[) will be established by induction on to, starting from to = fc. 
By ([4T5]| and flUTTJ we have, for any / G A k (T h ), 

tr/TTferfu = try^j^du + Py (du - n^du)} = dtr f P k ' 1 u + try (7 - Pf^^du. 

On the other hand, 

dtTf 7r£ _1 tt = dtxf Pf~ x u + dtry(7 - P^" 1 )^!^ 

= dti f P^ X u + try (7 - P^efor^M. 

By comparing the two expressions, and utilizing (|4.9[) , we obtain 

tr f (4du - drr^u) = tv f (I - P^)(4_ x d« - 

= try(7 - P^d^u - dnfclu) = 0, 

where the last identity is a consequence of Lemma T4. 61 So (|4.14[) holds for to = fc. 

Assume next that (|4.14[) holds for m replaced by to — 1. As we observed above, 
this implies that ir^ n _ 1 du = dn^Z^u. Furthermore, by (|4.12p it follows that the 
operators 7r^ _1 and tt^ satisfy (|4. 14[) for all / 6 A(T) with fc < dim/ < m — 1. 
Finally, for / e A m (T) we have by (O and (|4~TT|) that 

tr/ 7r^du = tr/[P*(du - ■K^ n _ l du) + ir^^du] 

= try dlPy*" 1 ^ - tt*" 1 ^) + Tr^n] = try d^ 1 "- 

This completes the proof. □ 
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5. Local bounds 

The purpose of this section is to derive local bounds for the projections tt^ 
constructed above. The main technique we will use is scaling, a standard technique 
in the analysis of finite element methods. The arguments below resemble parts of 
the discussion given in [1] Section 5.4], where scaling is used in a slightly different 
setting. 

^From the construction above, it follows that the operators ir k are local oper- 
ators. In fact, we observed in Section [3] that the operator ir%_i h = R k has the 
property that trfOTr k _ 1 h u only depends on u\ny As a consequence, (tt% _j h u)\r 
only depends on u restricted to 

U/eA fc (T)^/ C = A),r C Dk-i,r 

for T £ Tk and < k < n. Here we recall that the local domains D m T and 
D T = D n , T are defined by Therefore it follows by p3]) . (pL"TT]) . and the local 

properties of the operators P k and Ef, that the operator ir k has the property that 
(tt^u)\t only depends on u\d t for any T G Th, < k < n. Furthermore, for each 
h the operator 7r^ is a bounded operator in HA k (£l). Hence, for each h and each 
T G Th there is a constant c = c(h, T) such that 

(5.1) ||7r£u|| £aA » (T) < c{h,T) (\\u\\ L 2 Ak{DT) + \\du\\ L 2 Ak+HDT) ), ueHA k (D T ). 
Our goal in this section is to improve this result by establishing the uniform bound 

(5.2) ||7r£u|| iaAfc(T) <C7(|| u \\L 2 A k (D t ) + hT\\du\\ L 2 A k+i( DT )), u G HA k (D T ), 

for < k < n, where the constant C is independent of h and T . Since the operators 
Tr k commute with the exterior derivative, the estimate (|5.2[) will also imply that 

(5.3) ||d7r£u|| £3A » (T) <C\\du\\ L 2 Ak{DT) , ueHA k {D T ), 

for < k < n, with the same constant C as in (|5.2j) . Therefore, the estimate (|5.2j) 
will, in particular, imply the bounds given in Theorem 12.21 

The rest of this section will be used to prove the estimate (|5.2[) . For any fixed 
T £ Th, we introduce the scaling <3>t(^) = (x — Xo)/hr, where Xo is a vertex of T. 
We let T = $t(T) and Dt = &t(Dt) be the corresponding reference domains with 
size of order one. The restriction of the triangulation Th to Dt will be denoted 
Th{Dx), and Ti{Dt) the induced triangulation on Dt- In general we will use 
the hat notation to denote scaled versions of domains and local triangulations, for 
example / = $t(/), / € A(Th)- We note that the pullback, $ T maps HA k (D T ) 
to HA k (DT)- Furthermore, it follows from the definition of pullbacks that 

(5.4) WtAl^hd) =h T k+n/2 \H L 2 Ak{f)) , u€L 2 A k (D), 

where D C Dt and D = <&t{D)- We will obtain bounds for the operator 7r^, 
considered as a local operator mapping 7?A fe (Z>r) to HA k (T), by studying the 
operator ^ T ~ 1 n k i ^ T as an operator mapping HA k {r>T) to HA k {T). In fact, since 
the since the pullbacks commute with the exterior derivative, it follows from (|5.4[) 
that 

i k ii ll^*— 1 k n i — k+n/2 

lFhW|U=A"(T) = || $ T 7r fe U ll £ 2A"(T') /l r 
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(5.5) < pr^Tll ^"^(H^-^Wc?) + W^duh^nf)) 

< II^^V^^H (||u|| L 2 Afc(Dr) + h T \\du\\ L 2 Ak+HDT) ), 

where H^^rV^^H denotes the operator norm in C(HA k {D T ),L 2 A k (f)). Note 
that if we can show that this operator norm is uniformly bounded with respect to 
h and T E 7~h, then (|5.5j) will imply the desired bound (|5.2[) . The following result 
is the key tool for this verification. 

Lemma 5.1. The operator n k <fr T can be identified with the operator % k E 
C{H h k (Dt) , H k k (T)) obtained by constructing the operator Tr k with respect to the 
triangulation ThiDx) of Dt- 



Proof. We have to show that the operators n k and n k satisfies n k $ T = $> t tt . In 
fact, the proof just consist of checking that the pullback $ T commutes properly 
with the operators used to construct Tr k . A key property of the polynomial spaces 
VA k is that they are affine invariant. Therefore, in particular, we will have that the 
spaces VA k (Th(DT)) = $^'P A k (T}i(Dt)) ■ As a consequence of this we also obtain 
that the local projections Q k , defined with respect to the extended macroelements 
fiy, satisfies 

(5.6) $ T Q) = Qp* T , f E A k (%(D T )), 

with the obvious interpretation of Q k as the corresponding projections defined with 
respect to the domain f2* = $ T (^f)- A corresponding property holds for for the 
extension operators £ k , i.e., £ k (j) T = $ T £f> where £ k maps VoA k (f) to 7>f A k (Tf }h ). 
In particular, 

(5.7) Efvolf = $>* T £ k f <5>* T l sio\ f = <S>* T £ k f vo\j. 



Consider the operator S^ T , where S£ are the operators introduced in Section [3] 
above. By (|5 . T[) we have, for any u E HAr{pT), 



$tM A vol f2/ £f\/o\ f 



(5.8) 



S° h *5* = E { I 

f£A (T h (D T )) \ Jn f 
f£A (T h (D T )) V n f J 



/eA (T h (I>r)) 



In general, we define the operators S with respect to the reference domain Dt as 
outlined in Section[3] In particular, the weight functions z k are taken to be $ T z k . 
It follows essentially from (|5.6p . and an argument similar to one leading to 
that S k $ T = § T S k , and this further leads to 

(5.9) 4-i,h*T = Rh*T = *T%i = *T*k-l,h- 
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It is also straightforward to check that the local projections P k and the extension 
operators E k satisfy the corresponding properties P k <fr T = ^x^f anc ^ E k $ T = 
which implies that 

E) tTf Pf$T = ttf Pf. 

By combining this with the recursion ()4.11|) and ()5.9j) . we obtain the relation 
n m h^T = ^T^m f° r k <m< n. In particular, the desired relation 7r^$y = ^r^ k 
is obtained for m = n. □ 

We now have the following main result of this section. 

Theorem 5.2. The operators ir^ satisfy the bounds (|5.2j) and (|5.3|) . where the 
constant C is independent of h and T G 7/» • 

Proof. It follows from (|5.1|) that for each ft, and T, there is constant C(h,T) such 
that 

(5.10) ||7r fe U || L2Afe(f) <C(/i,T)H ffAfc(i j T)) ueHA k (D T ), 

where, as above, Tr k is obtained by constructing the operator ir k with respect to the 
triangulation Th{Dt) of Dt- However, due to the assumption of shape regularity 
of the family {7^}, it follows that the induced triangulations Th(Dt) varies over a 
compact set. Therefore, the constant C(h, T) is uniformly bounded with respect to 
h and T G 7fc. The desired estimate (|5.2p now follows from Lemma 15.11 combined 
with (|5.5[) and (|5.10|) . Finally, as we observed above, (|5.3[) follows from (|5.2[) and 
the fact that the projections nfc commutes with d. □ 

Finally, we observe that since the since shape regularity of the triangulation {Th} 
implies that the the covering {Dt} of 57 has a bounded overlap property, it follows 
from the bounds (|5.2|) and (|5.3|) that the global estimates 

lk^|| L 2 (0) < C (\\u\\ L 2 {n) + h\\du\\ L 2 m ) 

and 

\\dn^u\\ L 2 m < C(\\du\\ L 2 (n) , u G HA k (n), 
where C is independent of h, also holds. 
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